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ABSTRACT 

Foreground contamination is the fundamental hindrance to the cosmic microwave background (CMB) signals and 
its separation from it represents a fundamental question in Cosmology. One of the most popular algorithm used to 
disentangle foregrounds from the CMB signals is the "internal linear combination" method (ILC). In its original 
version, this technique is applied directly to the observed maps. In recent literature, however, it is suggested that in 
the harmonic (Fourier) domain it is possible to obtain better results since a separation can be attempted where the 
various Fourier frequencies are given different weights. This is seen as a useful characteristic in the case of noisy data. 
Here, we argue that the benefits of using such an approach are overestimated. Better results can be obtained if a classic 
procedure is adopted where data are filtered before the separation is carried out. 
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>- 1. Introduction 
l> 

f — . , The experimental progresses in the detection of cosniologi- 
fS| ■ cal emissions require a parallel development of data analysis 
' techniques in order to extract the maximum physical infor- 
' , mation from data. In particular, different emission mech- 
I ' anisms are characterized by markedly distinct underlying 
physical processes. Data analysis often requires the compo- 
, nent separation in order to study the individual character- 
IL" ' istics. To achieve such a goal, a link between the branch of 
^ \ signal processing science which characterizes and separates 
^ . different signals and astrophysics is yet well established, 
' and in many cases, modern signal processing techniques 
5^ [ have been imported and applied in an astrophysical con- 
text. 

In this work, we consider one of the most widely used 
approaches for the separation of different emissions, say 
the internal linear combination method (ILC), which was 
adopted, for instance, in the reduction of the data from the 
Wilkinson Microwave A nisotropy Probe (WM AP) satellite 
for CMB observations (jBennett et al. I [20031 ) . Among the 
separation techniques, ILC requires the smallest number of 
a priori assumptions. Here, the available data are assumed 
in the form of No maps, taken at different frequencies, con- 
taining Np pixels each. More precisely, if S'*^*)(p) provides 
the value of the pth pixel for a map obtained at channel 
" i " Q, our starting model is: 



where &i'\p), S^f'\p) and AA(*)(p) arc the contributions 
due to the CMB, the diffuse Galactic foreground and the 
experimental noise, respectively. Although not necessary 
for later arguments, it is assumed that all of these contri- 
butions are representable by means of stationary random 
fields. Moreover, without loss of generality, for ease of no- 
tation the random fields are supposed as the realization of 
zero-mean spatial processes. In the present work the contri- 
bution of non-diffuse components (e.g., due to SZ cluster, 
point-sources, . . . ) are not considered and they are assumed 
to have been removed through other methodologies. 

The idea behind ILC is simple. The main assumption is 
that model ((T|) can be written as 

S^^\p) = 6M + sl'\p)+Af^'\p), (2) 

i.e. the template of the CMB component is independent of 
the observing channel. A way to exploit this fact is to av- 
erage No images {5'^*^(p)}^i giving a specific weight Wi to 
each of the m so as to minimize t he impact of the foreground 
and noise (jBennett et al. II2003D . This means to look for a 
solution of type 

6e(p) = ^^/;,^«(p)- (3) 

i=l 

If the constraint X]i!=°i = 1 is imposed, Eq. ([3]) becomes 

e, (p) = s, ip) + [s^f^^ (p) + (p)] ■ (4) 

Now, from this equation it is clear that, for a given pixel 
"p" , the only variable terms are in the sum. Hence, under 



s^'Hp)^&^cHp) + s^;Hp)+^f^'Hp) (1) 



^ In the present work, p indexes pixels in the classic spatial 
domain. However, the same formalism applies if other domains 
are considered, for example, the Fourier one. 
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the assumption of independence of &cip) from Sf^\p) and 
J\f^''\p), the weights {wi} have to minimize the variance of 
&c{p), i-e. 



{wi} = argmin 



VAR [6c(p)] + VAR 



E 

.i=l 



(5) 



where VAR[s(p)] is the expected variance of s{p). 
If S^'^ denotes a row vector such as S^^' ~ 
[S'-'\l),S^'\2),...,S^''>{Np)] and the No x iVp matrix S 
is defined as 



then Eq. ^ becomes 
S 



©c + S'f +A/'. 



(6) 



(7) 



In this case, the weights are given by (jEriksen et al Il2004[ ) 

(8) 



w 



where Cs is the No x No cross-covariance matrix of the 
random processes that generate S, i.e. 



(9) 



and 1 = (1, 1, ... , 1)"^ is a column vector of all ones. Here, 
E[.] denotes the expectation operator. Hence, the ILC esti- 
mator takes the form 



with l^w = 1 and the scalar quantity a given by 



(10) 
(11) 

(12) 



In practical applications, matrix Cs is unknown and 
has to be estimated from the data. Typically, this is done 
by means of the estimator 



1 

N, 



-SS' 



(13) 



with 



/©c 

©1 

©2 



\&Nc/ 



(16) 



and 



Vi 



an ai2 



021 022 



ONoI ONa2 



02Nr_ 



(17) 



Here, matrix A is assumed to be of full rank. 



2. ILC in the spatial domain 

In a recent paper, IVio fc Andreani I (|2008l ) have stressed 
various problems concerning ILC that in literature have 
been underevaluated if not undetected. For example, in 
Eq. ([7]) the term Sf -t- A/" is often considered as a single noise 
component (e.g., see lEriksen et al71l2004t iHinshaw et al. I 
120071: iKim et al. Il2007l I2008D . In this way the problem is 
apparently simplified since it is reduced to the separation 
of two components only. No a priori information on this 
"global" noise is required. However, this approach can lead 
to wrong conclusions. For example, since all the compo- 
nents in the mixtures S are assumed to be zero-mean, from 
Eq. ([U one could conclude that 

E[©c|©c] = ©c + ti'^E[5f +A/'] = ©c, (18) 

i.e. the ILCjestimator is unbiased^ This is not correct: the 
claim that ©c is unbiased requires one to prove that 



E[©c|©c, S{] = ©c + w'^Sf + w'^E[J\f] = ©c 



(19) 



The reason is that 5'f is a fixed realization of a random 
process. There is no way to obtain another one. Even if 
observed many times (under the same experimental condi- 
tions) the foreground components (for instance the Galaxy) 
will always appear the same. Only the noise component 
A/" will change. This has im portant consequences. In fact, 
in IVio &: Andreani I (|2008[ ) it is proved that, even under 



model (|T5)l , ILC can provide satisfactory results only under 
rather restrictive conditions: 

1. The number of observing frequencies No has to be 
larger than the number of components N^, 



A common assumption in CMB observations is that 2. ©c has to be uncorrelated with 5'f; 

Sf'' is given by the linear mixture of the contribution „ „ , , , , . r ■ . /- 

' ^ , . , lAT / r r , 3. Data have to be noise-free, I.e. A/ = 

ot Nc physical processes {©jj^^i (e.g., tree-tree, dust re 

radiation, . . . ) 



] 1 



(14) 



with Oij constant coefHcients. In practice, it is assumed that 
for the jth physical process a template ©^ exists indepen- 
dent of the specific channel " i " . Although rather strong, 
it is not unrealistic to assume that this condition is satis- 
fied when small enough patches of the sky are considered. 
Inserting Eq. (fTil) into Eq. ([7]) one obtains 



The violation of any of these points has as consequence the 
introduction of a bias in the ILC solution that can be se- 
vere. In literature it appear s that only the second point ha s 
been well realized, (e.g. see iDelabrouille fc Cardoso1l2007[ ) . 
The explanation of the bias when No < Nc is technical 
and we remand to lVio fc Andreanil (|2008f) . When noise is 
present, under the condition of A/" uncorrelated with S, 
from model (fT5|) it is 



Cs = ACeA 



(20) 



S = A&+Af 



(15) 



^ The expression E[a|fe] indicates conditional expectation of a 
given 6. 
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with 



(21) 



Hence, a bias derives from the fact that Vlj^/ is a matrix 
with strictly positive entries and then E[Cs] Cs- 

In absence of "a priori inJormatior{\ the problems con- 
nected to the first two points above have no solution. The 
only possibility is to plan the experiments in such a way to 
avoid them (e.g. observations at high Galactic latitudes, a 
number of observing frequencies sufficiently large . . . ) . On 
the other hand, noise is an unavoidable question that, how- 
ever, can be handled with hope of satisfactory results. In 
this respect, son i e aut hors (e.g. see iTegmark et"ani2003t 
iKim et al. ll2007ll2008D suggest that an effective way is to 
use ILC in the harmonic (Fourier) domain. In practice, this 

means to apply ILC to S that is the matrix whose ith 
row conta ins the two-dimens ional Fourier transform of S*-*' . 
Following [Kim et al. I (|2007D . we will indicate this version 
of ILC as ^^harmonic internal linear combination^^ (HILC). 

3. ILC in the Fourier domain (HILC) 

The starting consideration is that all arguments presented 
in the previous sections hold independently of the fact that 
one is working in the ordinary spatial domain or in the 
Fourier domain. Indeed, it is not difficult to see that the 
weights w as given by Eq. ([5]) can be equivalently computed 
by means of 



w = 



(22) 



In this respect, it is sufficient to take into account that 



(23) 



with = ®Fm„, "fX"" the Kronecker product, Fjv the 
N X N the Fourier matrix that is a complex, symmetric and 
unitary matrix whose elements are given by 



1 



kl 



-27rt(fc-l)(i-l)/Af 



N 



-1. Hence 



(24) 



(25) 



since FF^ — p = j where F^ is the complex conju- 
gate transpose of Fat. 

It is well known that measurement noise {A/"^'-'}, typ- 
ically the realization of wide-band stochastic processes, 
tends to uniformly spread in the Fourier domain. On the 
other hand, the contribution of signals {S'*-*-'} is typically 
concentrated in correspondence to the lowest Fourier fre- 
quencies. For example, this is visible in Fig. [T] where the 
power spectra of three simulated mixtures are shown (see 
Fig. [3]). Here, non-astronomical signals have been used, 
but the same holds also for astronomical ones. This im- 
plies that, for the indices "p" corresponding the lowest fre- 
quencies, the contribution of noise to 5"^*^ (p) is quite small. 
Hence, the basic idea here is to partition the frequency do- 
main in a number Ns of subsets and to apply separately 
the ILC filter to each of them. At the end of this proce- 
dure an estimate &c = Ufc^i &c{Pk) is obtained where 



Pk ~ {p\p G k-th frequency subset}. Finally, &c can be 
recovered by the Fourier inversion of ©c- 

The effectiveness of such an approach seems supported 
by the simple example presented in Figs. [2E1 In particu- 
lar. Fig. [2] shows the original images &c, ©i, ©2 as well 
the mixing matrix A that, through model psp . are used 

to create the observed images {S^'^}l=i shown in Fig. [3 
As reference for later results, the bottom-right panel in the 
same figure shows the estimate &c obtained when ILC is 
applied to these images. For this case, the resulting root 
mean square (rms) of the residuals is about 0.05. The top 
panels in Fig. 3] show what happens when the images are 
added a white-noise with a signal-to-noise (SNR) set to 5. 
Here, both ILC and HILC are used. In the case of HILC, the 
frequency domain is partitioned in two regions, as shown in 
the bottom-left panel of Fig. [1] one containing the low fre- 
quencies and the other containing the high frequencies. The 
weights w are calculated independently for each of them. 
It is evident that HILC outperforms ILC. This conclusion 
is confirmed by the rms of the its residuals that is 0.13 
and 0.26, respectively. The same indication is provided by 
the bottom panels that show the result obtained when the 
weights computed for the noisy images are applied to the 
noise-free ones. This operation gives an idea of the bias in- 
troduced in the solution by the noise. In this case, the rms 
of the residuals become 0.08 and 0.18, respectively. 



4. Discussion and conclusions 

Reader should not be confused by the fact that at first 
sight HILC appears an approach potentially superior to 
ILC. It is indeed true that HILC performs a more effec- 
tive separation of the signal from the noise. However, if 
the noise affects most the high frequency part of the sig- 
nal, why then do not simply filter out this one? Indeed, 
the left panel of Fig. [5] shows the estimate &c obtained 
when, before the application of ILC, the images {S'''*') are 
filtered with the ideal Fourier circular low-pass filter □ that 
has the structure identical to that shown in Fig. [TJ As be- 
fore, the right panel of this figure shows the solution ob- 
tained when the resulting weights are applied to the noise- 
free version of the data. The rms of the residuals are 0.12 
and 0.08, respectively. The comparison with the values of 
0.13 and 0.08, previously determined for HILC, indicates 
an improvement of the quality of the result. The indication 
that comes out from this simple experiment is that filtering 
noise is a more effective operation than its separation from 
the signal of interest. It is true that the version of HILC 
here used is rather rough and m ore sophisticated ones are 
available (|Kim et al. II 2 00 7l [2008 ^ . However, the same holds 
also for the filtering operation that has been coupled with 
ILC. At this point a question arises: which is the benefit 
in using a non-standard approach as HILC with respect 
to a classic approach where filtering is coupled with ILC? 
This in not an academic question. Indeed, the use of non- 
standard tools implies the renunciation of a huge body of 
experience gained in years of application in practical prob- 



Here, the quantity SNR is defined as the ratio between the 
standard deviation of the signal with the standard deviation of 
the noise. 

* In the Fourier domain, an ideal filter H has the form H{p) = 
1 if p is a frequency of interest, otherwise. 
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lems also in disciplines different from Astronomy. In other 
words, the choice of non-standard tools is indicated only in 
situations of real and sensible improvements of the results. 
New techniques that do not fulfill this requirement should 
be introduced with care: they prevent the comparison with 
the results obtained in other works and may lead people to 
use not well tested methodologies ending up in not reliable 
results. 
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50 100 150 200 250 
Frequency 




50 100 150 200 250 
Frequency 



Power-spectrum N.3 (Log-scale) Partition of the frequency domain 




50 100 150 200 250 50 100 150 200 250 

Frequency Frequency 

Fig. 1. Top panels and bottom- left panel: power spectra (logarithmic scale) of the images S''-'-' shown in Fig. [3l The 
Figure shows that most of the power is concentrated at low frequency. Bottom-right panel: partition of the discrete 
Fourier domain that is used by HILC as discussed in Sect. [31 
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Original image N.1 (to be recovered) 



Original image N.2 
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Fig. 2. Top panels and bottom-left panel: original images &c, &i and ©2 - see Eq. (|16p - used in the experiment 
described in the text. The top-left panel shows the image to recover. Through Eq. ([TS]) these images are used to produce 

the mixtures S'*-'-' shown in Fig. [31 The bottom-right panel provides the mixing matrix A - see Eq. (I17p . 
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Observed N.I Observed N.2 




Fig. 3. Top panels and bottom-left panel: observed mixtures S'^^-', 5''^-' and S''^^ - see Eq. (|15l) - used in the experiment 
described in the text. The data shown here are noiseless. The right-bottom panel shows the solution obtained when ILC 
is applied to these images. 



8 



R. Vio, & P. Andreani: ILC method 



ILC solution (unfiltered images) HILC solution (all frequencies) 




ILC solution (noiseless case) HILC solution (noiseless case) 




Fig. 4. Top panels: solutions provided by ILC (left panel) and HILC (right panel) - see Sec. [3] when a white-noise, with a 
signal-to-noise ratio set to 5, is added to the images shown in Fig. [3] Bottom panels: Results obtained when the weights 
used to compute the solutions shown in the top panels are applied to the original (i.e. noise- free) images in Fig. This 
operation gives an idea of the bias introduced in the solution by the noise. 
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ILC solution (filtered images) ILC solution (noiseless case) 




Fig. 5. As in the left panels of Fig. [4] with the difference that before using ILC, the images have been filtered with the 
ideal circular low-pass filter that has the same structure as in the bottom-right panel of Fig. [TJ 



